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Considering gravitino dark matter scenarios, we study constraints on the reheating temperature 
of inflation. We present the gauge-invariant result for the thermally produced gravitino yield to 
leading order in the Standard Model gauge couplings. Within the framework of the constrained 
minimal supersymmetric Standard Model (CMSSM), we find a maximum reheating temperature of 
about 10 7 GeV taking into account bound-state effects on the primordial 6 Li abundance. We show 
that late-time entropy production can relax this constraint significantly. Only with a substantial en- 
tropy release after the decoupling of the lightest Standard Model superpartner, thermal leptogenesis 
remains a viable explanation of the cosmic baryon asymmetry within the CMSSM. 
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INTRODUCTION 

The observed flatness, isotropy, and homogeneity of 
the Universe suggest that its earliest moments were gov- 
erned by inflation [1, 2]. The inflationary expansion is 
followed by a phase in which the Universe is reheated. 
The reheating process repopulates the Universe and pro- 
vides the initial conditions for the subsequent radiation- 
dominated epoch. We refer to the reheating tempera- 
ture Tr as the initial temperature of this early radiation- 
dominated epoch of our Universe. 

The value of Tr is an important prediction of inflation 
models. While we do not have evidence for temperatures 
of the Universe higher than 0(1 MeV) (i.e., the temper- 
ature required by primordial nucleosynthesis), inflation 
models can point to Tr well above 10 10 GeV [2, 3]. 

In this Letter we consider supersymmetric (SUSY) ex- 
tensions of the Standard Model in which the gravitino G 
is the lightest supersymmetric particle (LSP) and stable 
because of R-parity conservation. The gravitino LSP is a 
well-motivated dark matter candidate. As the gauge field 
of local SUSY transformations and the spin-3/2 super- 
partner of the graviton, the gravitino is an unavoidable 
implication of SUSY theories including gravity [4]. Its 
interactions are suppressed by inverse powers of the (re- 
duced) Planck scale Mp = 2.4 x 10 18 GeV. Its mass mg 
results from spontaneous SUSY breaking and can range 
from the eV scale up to scales beyond the TeV region [5] . 

While any initial population of gravitinos must be di- 
luted away by the exponential expansion during infla- 
tion [6] , gravitinos are regenerated in scattering processes 
of particles that are in thermal equilibrium with the 
hot primordial plasma. The efficiency of this thermal 
production of gravitinos during the radiation-dominated 
epoch is sensitive to Tr [7-11]. Since the resulting grav- 
itino density £l~ p is bounded from above by the dark 
matter density Sld m , upper bounds on Tr can be de- 
rived [8, 12-15]. These bounds can be compared with 
predictions of the reheating temperature Tr from infla- 
tion models. Moreover, Tr is important for our under- 



standing of the cosmic baryon asymmetry. For exam- 
ple, successful standard thermal leptogenesis [16] requires 
Tr > 10 9 GeV [17]. 

We update the Tr limits using the full gauge-invariant 
result for the relic density of thermally produced graviti- 
nos, 0~ p , to leading order in the Standard Model gauge 
couplings [11]. This allows us to illustrate the depen- 
dence of the bounds on the gaugino-mass relation at the 
scale of grand unification Mqut — 2 x 10 16 GeV. 

We consider gravitino dark matter scenarios also in the 
framework of the constrained minimal supersymmetric 
Standard Model (CMSSM) in which the gaugino masses, 
the scalar masses, and the trilinear scalar interactions 
are assumed to take on the respective universal values 
m i/2i mo, an d Aq at Mqut- Taking into account grav- 
itinos from thermal production and from late decays of 
the lightest Standard Model superpartner, we provide 
new upper bounds on the reheating temperature in the 
(toi/2, mo) plane for various values of mg. Previous stud- 
ies of Tr constraints within the CMSSM used the result 
of [10] to explore the viability of Tr > 10 9 GeV [13, 14]. 
Our study presents also scans for Tr as low as 10 7 GeV 
based on the result of [11] which includes electroweak 
contributions to thermal gravitino production [18]. 

In the considered CMSSM scenarios with the gravitino 
LSP, the next-to-lightest SUSY particle (NLSP) is ei- 
ther the lightest neutralino Xi or the lighter stau t\} 
Because of the extremely weak interactions of the grav- 
itino, the NLSP typically has a long lifetime before it 
decays into the gravitino. If these decays occur dur- 
ing or after big-bang nucleosynthesis (BBN), the Stan- 
dard Model particles emitted in addition to the grav- 
itino can affect the abundance of the primordial light 
elements. Indeed, these BBN constraints disfavor the 
Xi NLSP for mg > 100 MeV [13, 14, 21]. For the 
slepton NLSP case, the BBN constraints associated with 



1 For simplicity, we consider Ao = in this work. For sizable |j4o|, 
also the lighter stop t\ can be the NLSP [19, 20]. 
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hadronic/elcctromagnetic energy injection have also been 
estimated and found to be much weaker but still signifi- 
cant in much of the parameter space [13-15, 21]. 

Only recently, it has been stressed that bound-state 
formation of long-lived negatively charged particles with 
the primordial nuclei can affect BBN [22-25]. With the 
charged long-lived stau NLSP, these bound-state effects 
also apply to the considered gravitino dark matter sce- 
narios. In particular, a significant enhancement of 6 Li 
production has been found to imply severe upper limits 
on the T\ NLSP abundance prior to decay [22] which 
strongly restricts the mass spectrum in the f\ NLSP 
case [26]. For generic parameter regions of the CMSSM, 
we show that this constraint disfavors Tr > 10 7 GeV and 
thereby successful thermal leptogenesis. 

Entropy production after decoupling of the NLSP and 
before BBN can weaken the BBN constraints signifi- 
cantly [27]. At the same time, the gravitino density is 
diluted which relaxes the bounds on Tr . We show explic- 
itly the effect of entropy production on the Tr bounds. 
Here we consider the cases of late-time entropy produc- 
tion before and after the decoupling of the NLSP. Indeed, 
a relaxation of the Tr bounds can render models of in- 
flation with T R > 10 7 GeV viable in CMSSM scenarios 
with gravitino dark matter. Since also a baryon asymme- 
try generated in the early Universe is diluted, the tem- 
perature required by thermal leptogenesis increases in a 
cosmological scenario with late-time entropy production. 
Still, we find that a sufficient amount of entropy produc- 
tion after NLSP decoupling and before BBN can revive 
successful thermal leptogenesis. 



THERMAL GRAVITINO PRODUCTION 

Gravitinos with rag > 1 GeV have decoupling temper- 
atures of T f G > 10 14 GeV, as will be shown below. We 
consider thermal gravitino production in the radiation- 
dominated epoch starting at Tr < T f G assuming that 
inflation has diluted away any initial gravitino popula- 
tion. 2 For T R < T f G , gravitinos are not in thermal equi- 
librium with the post-inflationary plasma. Accordingly, 
the evolution of the gravitino number density rig with 
cosmic time t is described by the following Boltzmann 



2 Taking a conservative point of view, we do not include gravitino 
production before the radiation-dominated epoch. However, in- 
fiaton decays, for example, can lead to a sizable yield of non- 
thermally produced gravitinos depending on the inflation model; 
cf. [28, 29] and references therein. 



TABLE I: The gauge couplings gi, the gaugino mass parame- 
ters Mi, and the constants a, ki, t/i, and (3^ associated with 
the gauge groups U(1) Y , SU(2) L , and SU(3) C . 
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equation [11] 

^+3Hn s = C 5 (1) 

where H denotes the Hubble parameter. The collision 
term Cg involves the gaugino mass parameters Mi, the 
gauge couplings gi, and the constants Cj and fcj associ- 
ated with the gauge groups U(1)y, SU(2)l, and SU(3) C 
as given in Table I. In expression (2) the temperature T 
provides the scale for the evaluation of Mj and gi. The 
given collision term is valid for temperatures sufficiently 
below the gravitino decoupling temperature, where grav- 
itino disappearance processes can be neglected. A pri- 
mordial plasma with the particle content of the minimal 
SUSY Standard Model (MSSM) in the high-temperature 
limit is used in the derivation of (2). 

The collision term (2) results from a consistent gauge- 
invariant finite-temperature calculation [11, 18] follow- 
ing the approach used in Rcf. [10]. Thus, in contrast to 
the previous estimates in [7, 8], the expression for Cg 
is independent of arbitrary cutoffs. Note that the field- 
theoretical methods of [30, 31] applied in its derivation 
require weak couplings, jj « 1, and thus high tempera- 
tures T > 10 6 GeV. 

Assuming conservation of entropy per comoving vol- 
ume, the Boltzmann equation (1) can be solved to good 
approximation analytically [10, 32]. At a temperature 
Tow -C Tr, the resulting gravitino yield from thermal 
production reads 

y tp (t \ = n j P ( r iow) Cg(r R ) 
G 1 lowj ~ s(T low ) ~ s(T R )H(T R ) 

- 2^ifc(rR)li + -3^-l 

i=l \ G / 

*Ksfe)(lo^v)' ,3) 

where the constants m are given in Table I. These 
constants are obtained with the Hubble parameter de- 
scribing the radiation-dominated epoch, iTad(T) = 
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^g*(T)7r 2 /90T 2 /Mp, the entropy density s(T) = 
27r 2 g*s(T) T 3 /45, and an effective number of relativis- 
tic degrees of freedom of <?*(Tr) = <7*s(Tr) = 228.75. 
We evaluate ^(Tr) and Mj(Tr) using the one-loop evo- 
lution described by the renormalization group equation 
in the MSSM: 



Mi(T) = 




-1/2 



M,(M GUT ) 



(4) 



(5) 



with the respective gauge coupling at the Z-boson mass, 
and the (3^ coefficients listed in Table I. 
Without late-time entropy production, the gravitino 
yield from thermal production at the present tempera- 
ture Tq is given by 



F| p (T )=^ TP (T low ) 



TP/ 

G 



(6) 



The resulting density parameter of thermally produced 
gravitinos is 



njh 2 = m d Y^(T ) s(T ) h 2 /p, 



TP , 



G 



"G^G 



(7) 



with the Hubble constant h in units of 100 km Mpc s 1 
and p c /[s(T a )h 2 } = 3.6 x 10~ 9 GeV. 

In Fig. 1 our result (3) for the thermally produced 
gravitino yield Y~ F (T\ ow ) is shown as a function of Tr 
for various values of rag (solid lines). The curves are 
obtained with m^n — 500 GeV for the case of univer- 
sal gaugino masses at Mgut: ^1,2,3 (Mqvt) — m i/2- 
The dotted lines show the corresponding results from the 
SU(3) C yield of Ref. [10] for M 3 = m 1 / 2l which was used 
to study Tr constraints on gravitino dark matter scenar- 
ios in Refs. [13-15]. We find that (3) exceeds the yield 
derived from [10] by about 50%; cf. [11]. The dashed 
(blue in the web version) horizontal line indicates the 
equilibrium yield 



G 



G 
s 



1.8 x 10 



-3 



(8) 



which is given by the equilibrium number density 
of a relativistic spin 1/2 Majorana fermion, n e ~ = 

3C(3)T 3 /(2tt 2 ). For T > if, fl »(T) = <?* S (T) = 230.75 
since the spin 1/2 components of the gravitino are in 
thermal equilibrium. In the region where the yield (3) ap- 
proaches the equilibrium value (8), gravitino disappear- 
ance processes should be taken into account. This would 
then lead to a smooth approach of the non-equilibrium 
yield to the equilibrium abundance. Without the back- 
reactions taken into account, the kink position indicates 
a lower bound for T f . 
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FIG. 1: The thermally produced gravitino yield (3) as a func- 
tion of Tr for m g = 10 MeV, 100 MeV, 1 GeV, 10 GeV, 
100 GeV, and 1 TeV (solid lines from left to right) and 
Mi i2 ,3(M G ut) = mi/2 = 500 GeV. The dotted lines show 
the corresponding yield obtained with the SU(3) C result for 
the collision term of Ref. [10]. The dashed (blue in the web 
version) horizontal line indicates the equilibrium yield of a 
relativistic spin 1/2 Majorana fermion. 



In the analytical expression (3) we refer to as the 
initial temperature of the radiation-dominated epoch. So 
far we have not considered the phase in which the coher- 
ent oscillations of the inflaton field dominate the en- 
ergy density of the Universe, where one usually defines 
Tr in terms of the decay width of the inflaton field <j). 
To account for the reheating phase, we numerically inte- 
grate (1) together with the Boltzmann equations for the 
energy densities of radiation and the inflaton field, 



dp r . 



fix I 



dt 



ad 



= r 



dt 



(9) 
(10) 



respectively; for details see Appendix F of Ref. [34] . 

With our result for the collision term (2), we find that 
the gravitino yield obtained numerically is in good agree- 
ment with the analytical expression (3) for 



Tr 



which satisfies T 



90 



-1 1/4 



3*(Tr)tt 2 



Ts M P 



1. 



(11) 



Towards smaller mg, T f G de- definition given by 



i-ffradC^k)- F° r an alternative Tr 

= £ -ffrad(TR), 



creases due to the increasing strength of the gravitino 
couplings. For example, for mg = 1 GeV (10 MeV), we 
find Tf > 10 14 GeV (10 10 GeV). 



90 



ff*(T R )7T 2 



1/4 



'r^Mp 



(12) 
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the associated numerically obtained gravitino yield is de- 
scribed by the analytical expression obtained after sub- 
stituting T R with y/£/1.8Tl in (3). 

While we focus on scenarios in which the gravitino is 
stable, the yield (3) is also crucial to extract cosmological 
constraints in scenarios with unstable gravitinos. Based 
on the result of [10] and taking into account thermal 
gravitino production during reheating, the following fit- 
ting formula was used to study constraints from decaying 
gravitinos in Refs. [33-35]: 



10 1 



Vf KM (T low ) ~ 1.9 x 1(T 12 



1 + 0.045 In 
1 - 0.028 In 



10 10 GeV 



10 10 GeV 

Tr 
10 10 GeV 



,(13) 



where Tr was defined via = 3iTad(TR). Compar- 
ing (13) with our result after the matching of the Tr 
definitions, we find that our result exceeds the mg- 
independent yield (13) by about 30% for m^, ;§> Mj(Tr). 
While the m^, dependence of Y~ p becomes negligible for 
decreasing Mj(TR)/mg, the yield (13) is used for 
as small as 100 GeV in Refs. [33-35]. As can be seen 
in Fig. 1, the actual yield for = 100 GeV is thereby 
underestimated by about an order of magnitude. Accord- 
ingly, the Tr bounds given in [33-35] are underestimated 
in the region < 1 TeV. 
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FIG. 2: Upper limits on the reheating temperature Tr. On 
the upper (lower) gray band, fi~ p for Mi, 2, 3 = rn 1 / 2 = 
500 GeV (2 FeV) at Mgut agrees with f^m- Fhe correspond- 
ing T R limits from the requirement ti-fh 2 < 0.126 shown by 
the dashed and dotted lines are obtained respectively with (3) 
for Mi/10 = M 2 /2 = M 3 = mi/2 at M G ut and with the re- 
sult of Ref. [10] for M 3 = m 1/2 at Mgut. 



CONSTRAINTS ON T R 

The reheating temperature Tr is limited from above 
in the case of a stable gravitino LSP since ~ p cannot 
exceed the dark matter density f^dm [8, 12-15]. In this 
paper, we use [36, 37] 

n&tf = o.io5l8;8& (14) 

as obtained from the measurements of the cosmic mi- 
crowave background (CMB) anisotropies by the Wilkin- 
son Microwave Anisotropy Probe (WMAP) satellite. 3 

In Fig. 2 we show the resulting upper limits on Tr 
as a function of mg. On the gray band, the thermally 
produced gravitino density (7) is within the nominal 3<r 
range (14). The upper (lower) gray band is obtained 
for M 12 ,3 = mi/2 at Mqut with m 1 / 2 = 500 GeV 
(2 TeV). The dashed lines show the corresponding con- 
straints for the exemplary non-universal scenario [39] 



This nominal 3cr range is derived assuming a restrictive six- 
parameter "vanilla" model. A larger range is possible — even with 
additional data from other cosmological probes — if the fit is per- 
formed in the context of a more general model [38] . 



Mi/10 = M 2 /2 = M 3 = m 1/2 at M GUT and from the 
requirement £l~ p ft. 2 < 0.126. Using the same require- 
ment and the result of [10] for M 3 = m^ 2 , we obtain the 
dotted lines. The electroweak contributions are particu- 
larly important for the considered case of non-universal 
gaugino masses at Mgut- For universal gaugino masses 
at Mqut, the Tr limits derived with the result of [10] 
provide already reasonable estimates. 4 

The Tr limits shown in Fig. 2 are conservative bounds 
that do only depend on mg and the values at Mgut- 
Once details of the SUSY model realized in Nature are 
known, one will be able to refine the limits by includ- 
ing contributions to fig from NLSP decays. In the next 
section, we will account for this non-thermal gravitino 
production in a systematic way within the framework of 
the CMSSM. 



4 Note that the dotted curve shown for m^ 2 = 500 GeV in Fig. 2 
is by about a factor of 4 more severe than the Tr limits shown 
in Figs. 5 and 6 of Ref. [14] in the region mg < 10 GeV in which 
fig P governs the limits. It seems to us that the gluino mass in 
(1.2) of Ref. [14] was accidentally evaluated at the scale fi = Tr 
rather than at the scale fj, ~ 100 GeV; see Sec. 5 of [10]. 
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CONSTRAINTS ON T R IN THE CMSSM 

In the CMSSM, one assumes universal soft SUSY 
breaking parameters at Mgut- The CMSSM yields phe- 
nomenologically acceptable spectra with only four pa- 
rameters and a sign: the gaugino mass parameter m 1 / 2l 
the scalar mass parameter uiq, the trilinear coupling Aq, 
the mixing angle tan/3 in the Higgs sector, and the sign 
of the higgsino mass parameter [i. 

Assuming Aq = for simplicity, the lightest Standard 
Model superpartner is either the lightest neutralino Xi or 
the lighter stau t\. Indeed, most CMSSM investigations 
assume that Xi is the LSP that provides dark matter; 
cf. [40] and references therein. The parameter region in 
which < m^o is usually not considered because of the 
severe upper limits on the abundance of stable charged 
particles [37] . In the gravitino LSP case, < m^o can 
be viable since the lightest Standard Model superpartner 
is unstable [13, 14, 25, 41]. 

With the gravitino LSP, the lightest Standard Model 
superpartner is the NLSP that decays into Standard 
Model particles and one gravitino LSP. For > 1 GeV, 
the NLSP decays after its decoupling from the thermal 
plasma. 5 Thus, the relic density of the associated non- 
thermally produced gravitinos reads [42] 

nl TP h 2 = m d Y NLSP (To) S (T )h 2 /p c (15) 

= S- ^NLSpft 2 , (16) 

"T-NLSP 

where tonlsp is the mass of the NLSP and Ynlsp(^o) and 
^nlsp'i are respectively the yield and the relic density 
that the NLSP would have today, if it had not decayed. 

In Fig. 3 the solid (black) and dotted (blue in the web 
version) lines show respectively contours of Inlsp(To) 
and tonlsp in the (rni/ 2 ,mo) plane for Aq = 0, ji > 0, 
tan/3 = 10 (left panel) and tan/3 = 30 (right panel). 
Above (below) the dashed line, rn^o < {m^ 1 < 

m^o). The medium gray and the light gray regions at 
small mi/2 are excluded respectively by the mass bounds 
m~± > 94 GeV and ran > 114.4 GeV from chargino and 
Higgs searches at LEP [37]. The leftmost dotted (blue 
in the web version) line indicates the LEP bound > 
81.9 GeV [37]. For tan/3 = 30, tachyonic sfermions oc- 
cur in the low-energy spectrum at points in the white 
corner labeled as "tachyonic." We employ the FORTRAN 
program SuSpect [43] to calculate the low-energy spec- 
trum of the superparticles and the Higgs bosons, where 
we use m t = 172.5 GeV for the top quark mass. Assum- 
ing standard cosmology, the yield Ynlsp(To) is obtained 



5 The NLSP freezeout temperature can be estimated from its 
mass: 7} NLSP < m NLSP /20 [12]. Thus, T R > T f NLSP for 
T R > 10 GeV which is considered in this Letter. 



from the ^nlsp^ 2 values provided by the computer pro- 
gram micrOMEGAs [44]. 

The contours shown in Fig. 3 are independent of 
and Tr. Therefore, they can be used to interpret the 
results shown in the figures below. Note the sensitivity of 
both (T ) and on tan (i. By going from tan (3 = 10 
to tan /3 = 30, (T ) decreases by about a factor of two 
at points that are not in the vicinity of the dashed line, 
i.e., that are outside of the ti-Xi coannihilation region. 
While m?! becomes smaller by increasing tan /3 to 30, the 
tan/3 dependence of m^o is negligible. 

Let us now explore the parameter space in which 

0.075 < nfh 2 + nf p h 2 < 0.126 . (17) 

Now, Tr and nig appear in addition to the traditional 
CMSSM parameters. We focus on mg > 1 GeV since 
the soft SUSY breaking parameters of the CMSSM are 
usually assumed to result from gravity-mediated SUSY 
breaking. However, we do not restrict our study to fixed 
relations between and the soft SUSY breaking pa- 
rameters such as the ones suggested, for example, by the 
Polonyi model. 

In Fig. 4 the light, medium, and dark shaded (green 
in the web version) bands show the (m 1/ / 2 ,mo) regions 
that satisfy the constraint (17) for Tr = 10 7 , 10 8 , and 
10 9 GeV, respectively, where tan/3 = 10, Aq = 0, \i > 0. 
The four panels are obtained for the choices (a) m^, — 
10 GeV, (b) m d = 100 GeV, (c) m d = 0.2 m , and 
(d) vtIq — rriQ. In the dark-gray region, the gravitino is 
not the LSP. The regions excluded by the chargino and 
Higgs mass bounds and the line indicating m^o = 
are identical to the ones shown in the left panel of Fig. 3. 
The dotted lines show contours of the NLSP lifetime. For 
the n NLSP, 

T fl cT- 1 (r 1 ^GT) = GLZ i__G (is) 

m\ \ m\ J 

as obtained in the limit m T — > 0. For the x? NLSP, we 
calculate t^o from the expressions given in Sec. IIC of 
Ref. [21]. 

The tnlsp contours in Fig. 4 illustrate that the NLSP 
decays during/after BBN. Successful BBN predictions 
therefore imply cosmological constraints on rog, mNLSP> 
and Ynlsp [13-15, 21]. Indeed, it has been found that 
the considered Xi NLSP region is completely disfavored 
for mg > 1 GeV by constraints from late electromag- 
netic and hadronic energy injection [13, 14, 21, 25]. In 
the Ti NLSP region, the constraints from electromag- 
netic and hadronic energy release are important but far 
less severe than in the x? NLSP case. Thus, much of 
the Ti NLSP region was believed to be cosmologically 
allowed [13-15, 21]. 

Recently, this picture has changed. It has been 
found that bound-state formation of long-lived negatively 
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FIG. 3: Contours of Ynlsp(To) (solid black lines) and ttinlsp (dotted blue lines) in the (mi/2,mo) plane for Aq = 0, n > 0, 
tan/3 = 10 (left panel) and tan/3 = 30 (right panel). Above (below) the dashed line, m-o < (m^ < rn~o). The medium 
gray and the light gray regions at small trii/j show the mass bounds m~± > 94 GeV and jtih > 114.4 GeV from chargino and 
Higgs searches at LEP [37]. 



charged T\ 's with primordial nuclei can catalyze the pro- 
duction of 6 Li significantly [22, 25]. Indeed, in most of 
the Ti NLSP parameter space, the associated bounds 
are much more severe than the ones from late energy 
injection. Only for < 10 3 s and > 40 GeV, 
the constraints from hadronic energy release can become 
more severe than the ones from catalyzed 6 Li produc- 
tion [25, 26]. We thus consider both the constraint from 
catalyzed 6 Li production derived in [22] and the one from 
late hadronic energy injection derived in [15]. 6 

For the constraint from bound-state effects on 6 Li pro- 
duction, we adopt the bounds given in Fig. 4 of Ref. [22] 
as Tfj -dependent upper limits on the yield of the neg- 
atively charged staus, Vnlsp/2. These bounds are ob- 
tained assuming a limiting primordial abundance of [45] 

( 6 Li/H) p < 2 x KT 11 . (19) 

The resulting constraint disfavors the t\ NLSP region to 
the left of the long-dash-dotted (red in the web version) 
line shown in Fig. 4. 

For the constraint from late hadronic energy injection, 
we use the upper limits on Inlsp that are given in Fig. 11 
of Ref. [15]. These limits are derived from a computation 



For details on the other BBN bounds and the additional CMB 
bounds, we refer the reader to the detailed investigations pre- 
sented in Refs. [13, 14, 21, 41, 46]. 



of the 4-body decay of the stau NLSP into the gravitino, 
the tau, and a quark-antiquark pair. 7 They are based on 
the severe and conservative upper bounds on the released 
hadronic energy (95% CL) obtained in [34] for observed 
values of the primordial D abundance of 

(WniOmean = (2.78±g;t|) x KT 5 (severe), 
(n D /n H )hi g h = (3.98t°;g7) x KT 5 (conservative). 

In Fig. 4 the associated constraints are shown by the 
short-dash-dotted (blue in the web version) lines. The D 
constraint disfavors the region between the corresponding 
lines in panel (b) and the region above the corresponding 
lines in panels (c) and (d) . In panel (a) the D constraint 
does not appear. 8 

Remarkably, one finds in each panel of Fig. 4 that 
the highest Tp, value allowed by the considered BBN 
constraints is about 10 7 GeV. The bands obtained for 
Tp > 10 8 GeV are located completely within the region 
disfavored by the 6 Li bound. In previous gravitino dark 



7 The 3-body estimate of the hadronic energy release given in 
Ref. [21] leads to overly restrictive limits, as shown in Ref. [15]. 

8 Additional constraints on hadronic energy release are imposed by 
the primordial abundances of 4 He, 3 He/D, 7 Li, and 6 Li/ 7 Li [25, 
34, 47-50]. However, in the region allowed by the 6 Li constraint 
from bound-state effects, i.e., Tf 1 < 10 3 s, the considered D 
constraint on hadronic energy release is the dominant one as can 
be seen in Figs. 38-41 of Ref. [34] and Figs. 6-8 of Ref. [50] . 
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(a) m 1/2 [GeV] (b) m 1/2 [GeV] 



m,Q = 0.2 m , tan /9 = 10, A = 0, ft > = m , tan /9 = 10, A = 0, /j, > 




100 200 500 1000 2000 5000 100 200 500 1000 2000 5000 

(c) m 1/2 [GeV] (d) m 1/2 [GeV] 



FIG. 4: The (m 1( i2,tno) planes for tan/3 = 10, Ao — 0, /i > 0, and the choices (a) mg = 10 GeV, (b) mg = 100 GeV, 
(c) m,Q = 0.2 mo, and (d) mg = mo. In each panel, the light, medium, and dark shaded (green in the web version) bands 
indicate the regions in which 0.075 < Q^h 2 < 0.126 for Tr = 10 7 , 10 8 , and 10 9 GeV, respectively. The medium gray and the 
light gray regions at small mi« are excluded respectively by chargino and Higgs searches at LEP. In the dark gray region, 
the gravitino is not the LSP. The dotted lines show contours of the NLSP lifetime. Below the dashed line, m^ < m^o. 
With the Ti NLSP, the region to the left of the long-dash-dotted (red in the web version) line is cosmologically disfavored by 
bound-state effects on the primordial 6 Li abundance [22]. The effects of late hadronic energy injection on the primordial D 
abundance [15] disfavor the ri NLSP region between the short-dash-dotted (blue in the web version) lines in panel (b) and the 
one above the corresponding lines in panels (c) and (d). The x? NLSP region above the dashed line, in which m^o < rrif 1 , 
is cosmologically disfavored by the effects of late electromagnetic/hadronic energy injection on the abundances of the light 
primordial elements [13, 14, 21, 25, 41]. 
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771^ = mo, tan/9 = 30, A = 0, ji > 




100 200 500 1000 2000 

m 1/2 [GeV] 



5000 



FIG. 5: As in Fig. 4 (d), for tan (5 = 30, A = 0, fi > 0, and 
mg = mo. 



matter studies within the CMSSM that did not take into 
account bound-state effects on the primordial 6 Li abun- 
dance, much higher temperatures of up to about 10 9 GeV 
were believed to be allowed [11, 13, 14]. 9 

The constraint Tr < 10 7 GeV remains if we consider 
larger values of tan (3. This is demonstrated in Fig. 5 for 
tan (3 = 30, Ao = 0, p > 0, and mg = jtiq. The shadings 
(colors in the web version) and line styles are identical to 
the ones in Fig. 4. 

Let us comment on the dependence of the considered 
BBN constraints on the assumed primordial abundances 
of D and 6 Li. As can be seen in Figs. 4 and 5, the con- 
straint from late hadronic energy release is quite sensitive 
on the assumed primordial D abundance. In contrast, 
even if we relax the restrictive 6 Li bound on Vnlsp/2 by 
two orders of magnitude, we still find Tr < 10 7 GeV. For 
example, the 6 Li constraint relaxed in this way would ap- 
pear in Fig. 4 (b) as an almost vertical line slightly above 
m 1/2 = 3 TeV. 

While the constraint Tr < 10 7 GeV is found for each 
of the considered mg relations, one cannot use the 6 Li 
bound to set bounds on without insights into mg. 
The 6 Li bound disappears for < 10 3 s [22] which is 
possible even for = O(100 GeV) provided is 
sufficiently small; see (18). However, the constraints on 
Tr become more severe towards small mg as is shown in 



Fig. 2. Thus, the constraint Tr < 10 7 GeV cannot be 
evaded by lowering mg provided Tr < T f G . 

An upper limit on Tr of 10 7 GeV can be problem- 
atic for inflation models and baryogenesis scenarios. This 
finding can thus be important for our understanding of 
the thermal history of the Universe. 



CONSTRAINTS ON T R WITH LATE-TIME 
ENTROPY PRODUCTION 

The constraints shown above are applicable for a stan- 
dard thermal history during the radiation-dominated 
epoch. However, it is possible that a substantial amount 
of entropy is released, for example, in out-of-equilibrium 
decays of a long-lived massive particle species X [2, 51]. 10 

If X lives sufficiently long, it might decay while its rest 
mass dominates the energy density of the Universe. The 
associated evolution of the entropy per comoving volume, 
S = sa 3 , is described by [2, 51] 



dS _ TxPxa 3 
dt ~ T 



27T 2 \ 1/3 

^5*) Ixpxa 4 ^ (20) 



together with the Boltzmann equation (10) for <f> = X 
and the Friedmann equation governing the evolution of 
the scale factor of the Universe a. Here Fx and px denote 
respectively the decay width and the energy density of X. 
Thus, the temperature after the decay can be expressed 
in terms of Tx, 



'aftc 



10 



#*(Tkfter)7r 2 



1/4 



(21) 



which satisfies Tx = 3if ra d(T a f tcr ). Indeed, primordial 
nucleosynthesis imposes a lower limit on this tempera- 
ture [55-58]: 



Tafter > 0.7-4 MeV 



(22) 



In Fig. 6 we show the evolution of S, a 3 px, and a 3 p ra d 
for two exemplary scenarios respecting (22). The scale 
factor a is normalized by a\ = a(10 GeV) = 1 GeV -1 and 
the temperature dependence of 5* is taken into account as 
determined in [59]. For p x (10 GeV) = 0.1p rad (10 GeV) 
and T a ftcr = 6 MeV, S increases by a factor of A = 100 as 
shown by the corresponding solid line. For px(10 GeV) = 
8/o ra d(10 GeV) and T a f tcr = 4.9 MeV, S increases by a 
factor of A = 10 4 as shown by the corresponding dotted 
(blue in the web version) line. 



Note that our bands for T R = 10 9 GeV differ from the ones 
shown in Refs. [13, 14]; see footnote 4. 



Gravitino dark matter scenarios with late-time entropy produc- 
tion have been considered previously for gauge-mediated SUSY 
breaking where T R > Tp [52-54] . 
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A = 100 
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FIG. 6: Evolution of S, a 3 px, and a 3 p ra d as a function of T 
for the normalization aj = a(10 GeV) = 1 GeV -1 . The solid 
lines are obtained for px(10 GeV) = 0.1p ra d(10 GeV) and 
Tatter = 6 MeV, the dotted (blue in the web version) lines for 
px(10 GeV) = 8p rad (10 GeV) and T aftcr = 4.9 MeV. 



We restrict our study to entropy production at late 
times, Tbcforc — T\ ow <C Tr, so that the thermal pro- 
duction of gravitinos is not affected. To work in a 
model independent way, we assume that the production 
of gravitinos and NLSPs in the entropy producing event, 
such as the direct production in decays of X, is negli- 
gible. 11 Moreover, in this section, we focus on scenar- 
ios in which the decoupling of the NLSP is not or at 
most marginally affected by entropy production, i.e., ei- 
ther T R » T after > t^lsp or prad ^ px for T > T rxLSP 

Thus, the thermally produced gravitino yield and — in the 
case of entropy production after NLSP decoupling — also 
the non-thermally produced gravitino yield are diluted: 



"J \-L after ) 



(23) 



In the case of late-time entropy production before the 
decoupling of the NLSP, we parameterize this by writing 



r| p (T ) = iy| p (r low ) 



(24) 



In this case, Vnlsp(To) and thereby fi| TP and the BBN 
constraints remain unaffected. 



11 The constraints discussed below shall therefore be considered as 
conservative bounds. For studies of gravitino production dur- 
ing an entropy producing event, we refer to [60] and references 
therein. 



Conversely, in the case of late-time entropy production 
after the decoupling of the NLSP (and before BBN) both, 
Y| P (T ) and Y N lsp(?o), are reduced: 



Y? P (T ) 



1 
1 



TP 



low ) 



NLSP 



(To) — ^iNLSp(^low) 



(25) 



become smaller and the 



Accordingly, Sl~ p and f2~ TP 
BBN constraints can be relaxed. 

In Fig. 7 we show how late-time entropy production be- 
fore (left) and after (right) NLSP decoupling affects the 
6 Li constraint and the region in which 0.075 < fi^h 2 < 
0.126 for Tr = 10 9 GeV. The (toi/ 2 , m ) planes are con- 
sidered for tan/3 = 10, A = 0, fi > 0, = 100 GeV 
(upper panels) and = too (lower panels). The dark 
shaded (green in the web version) region is obtained with- 
out late time entropy production, 5 = A = 1. The 
medium and light shaded (green in the web version) 
bands are obtained with a dilution of ~ p (^ P + ^^ TP ) 
by 5 = 10 (A = 10) and 6 = 100 (A = 100), respectively. 
The dot-dashed (red in the web version) line illustrates 
that the 6 Li bound is independent of S, as shown in the 
panels on the left-hand side, and becomes weaker (i.e., 
moves to the left) with increasing A, as shown in the 
panels on the right-hand side. Other curves and regions 
are identical to the ones in the corresponding panels of 
Fig. 4. Note that we do not show the D constraint on 
late hadronic energy injection since it is not sensitive to 
S and vanishes already for A = 10; an exception is the 
severe D constraint which still appears for A = 10 in 
panel (d). BBN constraints on Xi NLSP scenarios with 
entropy production after NLSP decoupling will be stud- 
ied elsewhere. 

Comparing panels (b) and (d) of Fig. 4 with panels (a) 
and (c) in Fig. 7, we find that a dilution factor of S = 10 
(100) relaxes the Tr bound by a factor of 10 (100). Since 
the BBN constraints are unaffected by 5, the cosmolog- 
ically disfavored range of NLSP masses cannot be re- 
laxed. With the dilution after NLSP decoupling, the re- 
laxation of the Tr constraints is more pronounced. Here 
also the cosmologically disfavored range of NLSP masses 
can be relaxed [27]. However, as can be seen in panels (b) 
and (d) of Fig. 7, the 6 Li bound is persistent. With a di- 
lution factor of A = 100, large regions of the (toi/2,too) 
plane remain cosmologically disfavored. For A > 10 4 , 
however, the 6 Li bound can be evaded as will be shown 
explicitly below. 

Figure 7 shows that inflation models predicting, for 
example, T R = 10 9 GeV become allowed in the CMSSM 
with gravitino dark matter for S = A s» 100. Here it 
is not necessary to have late-time entropy production in 
the somewhat narrow window between NLSP decoupling 
and BBN. This is different for the viability of thermal 
leptogenesis in the considered scenarios (T f G > Tr) and 
for collider prospects as discussed below. 
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= 100 GeV, tan /3 = 10, A = 0, ^ > = 100 GeV, tan /3 = 10, A = 0, ^ > 
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FIG. 7: The effect of fate-time entropy production before (feft) and after (right) NLSP decoupling on regions in which 0.075 < 
fl Q h 2 < 0.126 for T R = 10 9 GeV. The (m 1/2 ,m ) plane is shown for tan/3 = 10, Ao = 0, fx > 0, m & = 100 GeV (upper 
panels) and mg = mo (lower panels) . The dark shaded (green in the web version) region is obtained without late-time entropy 
production 5 = A = 1. The medium and light shaded (green in the web version) bands are obtained with a dilution of f2~ p 
(ng p +0| TP ) by 8 = 10 (A = 10) and S = 100 (A = 100), respectively. The fi NLSP region to the right of the dot-dashed (red 
in the web version) line is cosmologically disfavored by the primordial 6 Li abundance. Other curves and regions are identical 
to the ones in the corresponding panels of Fig. 4. The severe D constraint for A = 10 appears only in panel (d) . 



THERMAL LEPTOGENESIS IN THE CMSSM tory strongly disfavors thermal leptogenesis. However, 
WITH GRAVITINO DARK MATTER if entropy is released after NLSP decoupling, a dilution 

factor of A ~ 10 4 can render thermal leptogenesis viable 
The constraint T R < 10 7 GeV obtained in the consid- for T R ~ 10 13 GeV. 
ered CMSSM scenarios for a standard cosmological his- 
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Standard thermal leptogenesis usually requires Tr > 
10 9 GeV [17]. However, late-time entropy production 
dilutes the baryon asymmetry which is generated well 
before NLSP decoupling, 

^(Tafter) = i ^(T bc foro) ■ (26) 

Therefore, the baryon asymmetry before entropy pro- 
duction must be larger by a factor of A in order to 
compensate for the dilution. For A ~ 10 4 , this can 
be achieved in the case of hierarchical neutrinos for 
Mri ~ T R ~ 10 13 GeV, as can be seen in Fig. 7 (a) of 
Ref. [61] and in Fig. 2 of Ref. [62]. Here M m is the mass 
of the lightest among the heavy right-handed Majorana 
neutrinos. 

In Fig. 6 the dotted (blue in the web version) lines 
show a scenario in which a dilution factor of A = 10 4 
is generated in the out-of-equilibrium decay of a heavy 
particle X. Because of px(10 GeV) = 8p ra d(10 GeV), the 
Hubble rate can be enhanced already during the decou- 
pling phase of the NLSP, which leads to an increase of 
T f NLSP and V N LSp(J 1 f NLSP ). In the results shown below, 
we account for this by using a modified version of the 
micrOMEGAs code. 12 After entropy production, the net 
effect is still a significant reduction of Fnlsp(To). For 
the same initial conditions, A = 2 x 10 4 — and thereby an 
additional reduction of Inlsp(^o) by a factor of two — 
can be achieved by lowering T a ftcr from 4.9 MeV down to 
2.5 MeV. 

We consider these two scenarios for tan /3 = 30, A = 0, 
/j > 0, and rrig = too, in Fig. 8. Here the shaded (green 
in the web version) bands indicate the region in which 
0.075 < ft d h 2 < 0.126 for T R = 10 13 GeV and A = 10 4 
(dark) and 2 x 10 4 (medium). In addition, the corre- 
sponding evolution of the 6 Li bound is shown by the dot- 
dashed (red in the web version) lines. For A = 10 4 , the 
regions below the associated two rightmost curves and 
to the right of the associated leftmost curve are allowed. 
For A = 2 x 10 4 , the cosmologically allowed region is the 
Ti NLSP region below the line labeled accordingly. The 
gray regions are identical to the ones in Fig. 5. 

We find that the 6 Li bound cannot be evaded for the 
tan (3 — 10 scenarios even for A = 2 x 10 4 since Inlsp(To) 
becomes larger. However, the 6 Li bound given in Fig. 4 
of Ref. [22] depends linearly 13 on the assumed limiting 
primordial abundance (19) that is subject to uncertain- 
ties; cf. Ref. [50]. Accordingly, for a limiting abundance 
that is a factor of two above the value given in (19), one 
obtains the 6 Li bound labeled with A = 2 x 10 4 in Fig. 8 
for the scenario with tan/3 = 30 and A = 10 4 . 

Scenarios with successful thermal leptogenesis in the ti 
NLSP region are located preferably on the dark-shaded 



The Ynlsp contours shown in Fig. 3 do not apply in this section. 
We thank M. Pospelov for bringing this point to our attention. 



nig = mo, tan /3 = 30, A = 0, fi > 




100 200 500 1000 2000 5000 

m 1/2 [GeV] 

FIG. 8: The effect of entropy production after NLSP decou- 
pling for T R = 10 13 GeV and A > 10 4 in the (m 1/2 ,m ) 
plane for tan/9 = 30, Ao = 0, /i > 0, and — mo. The 
shaded (green in the web version) bands show the region in 
which 0.075 < tt d h 2 < 0.126 for A = 10 4 (dark) and 2 x 10 4 
(medium). The dot-dashed (red in the web version) lines il- 
lustrate the corresponding evolution of the 6 Li bound. For 
A = 10 4 , the regions below the associated two rightmost 
curves and to the right of the associated leftmost curve are 
allowed. For A = 2 x 10 4 , the region below the line labeled 
accordingly is cosmologically allowed. 



(dark green in the web version) band and in the white 
corner to its left, in which even slightly higher values of 
T R are possible for A = 10 4 . For T R = 10 13 GeV and 
A ^> 10 4 , the generated baryon asymmetry is diluted too 
strongly in order to explain the observed baryon asym- 
metry. 

As can be seen in Fig. 8, the t\ NLSP region with 
500 GeV < m 1/2 < 700 GeV, where m ?1 < 200 GeV 
(cf. Fig. 3), is no longer disfavored by the 6 Li bound pro- 
vided A > 10 4 . Such scenarios are particularly promising 
since the long-lived T\ NLSP could provide striking signa- 
tures of gravitino dark matter at future colliders [63-66] . 



CONCLUSION 

Using the full gauge-invariant result for il~ p to lead- 
ing order in the Standard Model gauge couplings [11], 
we have studied bounds on Xr from the constraint f2g < 
fidm- Our results take into account the dependence of 
O ~ p on the masses of the gauginos associated with the 
Standard Model gauge group SU(3) c xSU(2) L xU(l) Y . 
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This has allowed us to explore the dependence of the 
Tr bounds on the gaugino-mass relation at the scale of 
grand unification Mqut- 

Within the CMSSM, we have explored gravitino dark 
matter scenarios and the associated Tr bounds for mg > 
1 GeV and for temperatures as low as 10 7 GeV. Taking 
into account the restrictive constraint from bound-state 
effects of long-lived negatively charged staus on the pri- 
mordial 6 Li abundance [22], we find that T R < 10 7 GeV 
is the highest cosmologically viable temperature of the 
radiation-dominated epoch in case of a standard thermal 
history of the Universe. This imposes a serious constraint 
on model building for inflation. Moreover, thermal lepto- 
genesis seems to be strongly disfavored in the considered 
regions of the CMSSM parameter space. 

With late-time entropy release, the obtained limit 
Tr < 10 7 GeV can be relaxed. For example, the dilution 
of the thermally produced gravitino yield by a factor of 
10 relaxes the Xr bound by about one order of magni- 
tude in regions where f2~ p dominates fig. In the case of 
entropy production after NLSP decoupling, the yield of 
the NLSP prior to its decay, VnlsP: is reduced so that 
the BBN constraints can be weakened. Although the 6 Li 
bound is persistent, we find that it disappears provided 
Vnlsp is diluted by a factor of A > 10 4 . 

We have discussed the viability of thermal leptogen- 
esis in a cosmological scenario with entropy production 
after NLSP decoupling. We find that successful ther- 
mal leptogenesis can be revived in generic regions of the 
CMSSM parameters space for M m ~ T R ~ 10 13 GeV 
and A > 10 4 , where Mri is the mass of the lightest 
among the heavy right-handed Majorana neutrinos. 

Remarkably, for a dilution factor of A > 10 4 , the t\ 
NLSP region with < 200 GeV reopens as a cosmolog- 
ically allowed region in the CMSSM with the gravitino 
LSP. A long-lived t\ in this mass range could provide 
striking signatures of gravitino dark matter at future col- 
liders [63-66]. 

We are grateful to S. Blanchct, F. Hahn-Woernlc, 
M. Plumachcr, M. Pospelov, G. Raffelt, and 
Y.Y.Y. Wong for valuable discussions. 
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